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In order to explore various aspects of stellar evolution, supernovae, gamma ray bursts 
and nucleosynthesis, we have developed a new efficient stellar evolution code. In this paper 
we describe this new code and compare the results with the ones calculated by the previous 
code. Specifically we focus on the progenitor evolution of lower end of the Fe-core collapse 
supernovae, and mass distribution of remnant neutron stars. We describe how different 
assumptions will lead different neutron star mass distribution. We also review recent works 
of our research group. 

§1. Introduction 

r~| ! Massive stars end their life as supernovae leaving neutron stars behind, or by 

forming blackholes without explosions if they are not rotating. The critical initial 
stellar mass for the blackhole formation is usually considered to be about 20 to 25 
Mq but it is uncertain.'!'' This is because explosion mechanism for supernovae is 
not yet known definitely]^ ''^'^''^''^''^ and also there are still uncertainties in the 
progenitor models. 

Lower end for the neutron star forming supernova is also still uncertain (e.g., 
Ref. IHj)). Here, the uncertainties in the stellar evolution theory may be even larger. 
It is well known that the stars above around 10 Mq stars form an Fe core in the end 
of their evolution. The formation of Fe core is relatively simple for a M > \2>Mq 
star (review in Ref. [9])). In such a star the hottest region is the center mostly all the 
time through the evolution. Therefore, heaviest element is synthesized around the 
center, forming Fe-core in the end. 

On the other hand the evolution of M < ISMq stars are more complicated. 
Less massive stars tend to have temperature inversion between the center and outside 
regions. This phenomenon is well known for intermediate mass stars with M < 8Mq 
after carbon burning, for which the carbon burning starts at off center (e.g., Ref. llOp ). 
For a M < IIMq star, off center O- and/or Si-burning occurs. The off-center 
burning front propagates inward through complicated burning stages, forming an 
Fe core eventually. This off center burning becomes sometimes very violent and it 
could cause some mass ejectiorpJ-Eli (^ggg a,lso Section 3.1 in this paper). However, 
such calculations have not been done recently with updated input physics, so it is 
currently not clear if such mass ejection really occurs. 

For less massive stars which do not ignite Ne even at off center, a cool degenerate 
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0-Ne core is first formed. Tliis 0-Ne core grows in mass gradually when the star is 
in a super AGB phase. Depending on the mass loss rate this 0-Ne core reaches the 
critical mass for core-collapse owing to electron capture.^^® It is considered that 
such a star explodes by the neutrino energy transport mechanism.^ 

Calculations for the progenitors of such low mass core-collapse supernovae with 
updated input physics are interesting and important. However, their evolution can 
be quite different if the initial mass is only slightly different, say by 0.01-0.1 Mq 
(e.g., Ref. i8j); see also Sec. 3.1). To understand the whole story in this mass rage, 
therefore, we need to calculate stellar evolution in very fine mass grids. This will be 
quite time consuming and thus currently well-used progenitor models in the literature 
(e.g., Refs.[I2]),[I5I),[IS]),[IZD,[ISI),[I2])) do not fully deal with this mass range. 

To tackle this problem, we have developed a stellar evolution code for efficient 
computation. In this paper we describe this new code, the Yoshida-Umeda (YU) 
code^ and compare the results with the ones calculated by one of the author, 
H.U., using the Umeda-Nomoto (UN) code^^'^^ We also briefly review some 
of the recent works using the YU code and other works in our research group. 

This paper is organized as follows. In §2, we describe the new stellar evolution 
code and differences with previous codes. In §3, massive star evolutions calculated 
with this code are compared with the ones with a previous code. In this section we 
also describe the progenitor evolution of lower end of Fe-core collapse SNe in some 
detail. §4 describes nucleosynthetic aspects including some reviews of our recent 
works. In §5 remnant neutron star masses are given as a function of progenitor 
mass. We describe how the different assumption will lead different neutron star 
mass distribution. §6 reviews other recent works of our research group and §7 gives 
discussions and Future Prospects. 

§2. Stellar evolution code 

2.1. Umeda-Nomoto (UN) code 

Before describing the YU code, we briefly describe the UN code because we will 
compare the results of these codes. This code is mostly based on Nomoto-Hashimoto 
(NH) code^^' and Saio-Nomoto-Kato (SNK) code.'^ In the UN code, input physics 
such as equation of state (EOS) are same as the NH code except updates of radiative 
opacity, neutrino emissivity and electron capture rates. 

There are several differences in the NH and SNK codes: 

i) The NH code is a He-star code, which means that it cannot solve stellar 
atmosphere and hydrogen burning phase. Thus it cannot provide stellar radius 
correctly. On the other hand, the SNK code solves atmosphere. 

ii) Since the SNK code is not designed for calculating later evolutionary stages 
of massive stars, it does not deal with nuclear burning after carbon burning stages. 
Also because of that, the SNK code does not include the acceleration terrrPS and 
omit the inertial term in the equation of motion, while the NH code can include 
it. This acceleration term becomes important for constructing supernova progenitor 
models just before iron-core collapse. 
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iii) The treatment of convection is different. Both codes adopt Schwarzschild 
criterion for convection but the way of mixing is different. The NH code assumes 
instantaneous mixing in convective regions, while in the SNK code matter is mixed 
diffusively using the formalism of Ref. 1251) taking into account of the semi-convection 
effects. The energy transfer in a convective region is also different. In the NH code 
time-dependent mixing length theorj^'^^S can be included for convective energy 
transfer, though this effect is not so important for the massive star evolution forming 
an Fe-core. 

In the UN code atmosphere is calculated as the SNK code, the acceleration 
term is included and the formalism of Ref. [25]) is adopted for convective mixing. 
Convective energy transfer is treated as same as the NH code. 

The UN code differs from the NH and SNK codes in the calculations of nucleosyn- 
thesis and nuclear energy generation. The NH code assumes quasi nuclear statistical 
equilibrium during the Si burning while the UN code solves full nuclear reaction net- 
works below logr(K) < 9.6. Above that temperature both codes assume nuclear 
statistical equilibrium (NSE). In order to calculate nuclear energy generation rates, 
in the UN code nuclear reaction networks are solved simultaneously with Henyey 
relaxation, while in the NH and SNK codes abundance is fixed during Henyey relax- 
ation. In this sense the UN code solves abundance implicitly, while the NH and SNK 
codes solve explicitly. Solving abundance implicitly is the best way to obtain con- 
sistency in the energy generation rates and the abundance evolution. However, this 
has a disadvantage in efficient calculations because solving large reaction networks 
involves time consuming matrix inversion calculations. 

2.2. Yoshida & Umeda (YU) code 

Basic structure of this code is based on the SNK code.'^ As mentioned above, 
although the UN code has been successfully used for the progenitor calculations, it 
has an disadvantage in the calculation time. To finish a calculation from ZAMS to 
Fe-core collapse it typically takes few months. Therefore this code is not suitable 
for a large parameter search. Hence, we have developed a new more efficient code 
(YU code). Main difference of the YU code from the UN code is the calculations of 
nucleosynthesis and energy generation. 

The YU code solves a full nuclear reaction network from hydrogen burning up 
to log^o Tq ~ 10.0. The nuclear reaction network consists of ~ 300 species of nuclei 
from n, p to Br. NSE is not assumed in calculations. On the other hand, the 
nucleosynthesis is solved before Henyey relaxation in the YU code similar to the 
SNK code. Using the abundance of the next step, the energy generation is calculated 
during the Henyey relaxation. Therefore, smaller time step is required. 

After carbon burning stage, the time interval is typically determined to satisfy 
the conditions Z\ log^o logio ^ ^ 0.001 and log^g p/ logio p ^ 0.003 for each 
mass coordinate. This treatment requires more calculation steps but still saves the 
calculation time especially in the late phase of the evolution. With this code, we can 
calculate one model typically in 1 — 2 weeks for massive star evolutions. 

This code at present has not implemented the acceleration term yet. We will 
include that term in near future. 
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Fig. 1. Evolution tracks of M = 12, 15, 20, and 25AIq stars in terms of central temperature and 
density. Mass fraction distributions of a 15Mq star denoted at the points a — e are shown in 
Fig-H 

2.3. Mass loss rates 

In the UN and YU codes, same mass loss rates are used for the solar metallicity 
before the Wolf-Rayet star phases shown in this paper. In OB stars, where the surface 
temperature is larger than 1.2 x 10^ K and the surface H mass fraction is larger than 
0.4, the mass loss rate is adopted from Ref. [28]). In yellow supergiants and red- 
giant branches, the mass loss rate in Ref. [29j) is used. In the UN code, metallicity 
dependent factor of is multiplied to the rate for metal poor stars as 

Ref. [30|1 . The YU code adopts the metallicity dependence in a different manner. 
The metallicity dependence in the main-sequence stage is taken from Ref. [28]). In 
the YU code presented in this paper the case A mass loss rate in the Ref. [20j) is 
used. In yellow supergiants and red-giant branches, the metallicity dependent factor 
(Z/0.02)'^'^^ is multiplied to the mass loss rate. The power index is the same as that 
of B supergiants in Ref. I28]) . See Ref. I20]) for the detail. 

§3. Progenitor Evolution 

Basic properties of massive star evolution have been discussed in many papers 
(e.g., Refs. [T2]) . [T6]) .[T7 )) .[T8 ]) .[T9 ]) ). therefore we do not repeat detail here, except 
M ~ 10 — 13Mq models in the next subsection. We briefly describe the evolution in 
advanced stages of a M = 15Mq star as an example. Fig. [T]shows the evolution tracks 
of M = 12, 15, 20, and 25Mq stars in terms of the central density and temperature. 
We also show the mass fraction distributions of the 15 Mq star at five stages during 
the evolution in Fig. [21 Fig. ^a) indicates the mass fraction distribution after the 
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Fig. 2. Mass fraction distributions of a i\f = ISAf© star after the He burning (a), C burning (fe), 
Ne burning (c), O burning (d), and Si burning (e). The final mass is Mi — 13.6A4q. "Si" and 
"Fe" indicate the sum of the mass fractions of the elements Si-Sc and Ti-Br. 



core helium burning. There is a CO core of 2.2 Mq in the central region. The 
He layer of 2 Mq and the H-rich envelope surround the CO core. The He layer 
consists of ashes of the hydrogen burning. After the core helium burning, the shell 
helium burning produces ^^C in the He layer. Then, the carbon ignites at the center 
when the central temperature becomes log^g Tc(K) ~ 8.8. The core carbon burning 
converts ""^^C into ^"^Ne and ^'^Na to form an 0/Ne core (Fig. E^^)). Some Mg-Si are 
also produced in the core. 

The shell carbon burning is followed by the core carbon burning. It converts the 
0/C layer into an 0/Ne layer. However, ^^C of 0.09 by mass fraction remains even af- 
ter shell carbon burning. When the central temperature becomes logj^g Tc(K) ~ 9.2, 
the core neon burning starts and an 0/Si core forms (Fig. [2Jc)). After the Ne burn- 
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ing, oxygen ignites at the center with the central temperature of log^g ^c(K) ~ 9.3. 
The core oxygen burning produces a ~ IMq Si core (Fig. E^d)). The main com- 
ponents of the Si core is ^^Si, '^■^S, and '^^S and the electron fraction decreases by 
electron captures. The Si core extends through the following shell oxygen burnings 
and the central temperature increases by the contraction. When the central temper- 
ature is logio7c(K) ~ 9.55, the Si burning starts and an Fe core forms (Fig. EJe)). 
The composition in the O layer also gradually changes through shell burnings. After 
the Fe core formation, the Fe core extends with shell Si burning and finally collapses 
to explode as a supernova. 

Here we stress that one of the most important factors in the massive star evolu- 
tion is the carbon abundance after core helium burning, Xc(^^C). In general smaller 
carbon abundance leads a larger iron core at core collapse. This is because when the 
carbon abundance is smaller, shell carbon burning is weaker after helium burning. 
Then the core is less supported by the convective carbon shell burning and the time 
between the central helium burning to the oxygen burning is shorter. As a result, 
for smaller carbon abundance a larger core with higher entropy is formed because 
neutrino cooling is less effective .'^''^ 

The carbon abundance also affects on nucleosynthesis substantially. When the 
abundance is too low, carbon burning products such as Ne, Na, Mg, and Al are 
underproduced compared with the solar abundance pattern, while they are overpro- 
duced when the abundance is too large.l^D They also found that the abundances of 
S, Ar, and Ca are anti-correlates with Ne, Na, Mg and Al. 

Despite of these importance the carbon abundance is very uncertain because 
it sensitively depends on the ^^C(a,7)^^0 rate and the treatment of convection. 
Ref. I3T]) discussed that due to the nucleosynthetic constraint described above, in 
their model the ^■^C(a,7)^^0 rate needs to be 1.7 ± 0.5 times the Caughlan and 
Fowler (1988) (CF88 hereafter) rate.l^^^ In the NH88 modeP* the rate of Ref. [MD 
was taken which corresponds to 2.3 - 2.4 times the CF88 rate. In NH88, Xc(^^C) = 
(0.25, 0.22, 0.19) for the helium star model for Mne= (3.3, 6, 8) Mq, which roughly 
corresponds to the ZAMS M=(13, 20, 25) Mq, respectively. In the UN model the 
"'^^0(0,7)^^0 rate was chosen to be 1.3 times larger than the CF88 rate so that the 
abundance pattern of EMP and VMP stars are reproduced.'22''El 

Although there have been improvements in the estimation of the ^^C(a;, 7)^^0 
uncertainty in the rate is still about factor of two. Therefore, the nucle- 
osynthetic method is still the best way to constrain the carbon abundance. 

We show in Table I the central carbon abundance after the core helium burning 
for UN and YU models. In the YU model ^^C(a,7)^^0 rate is chosen to be 1.5 
times as large as the CF88 rate so that the carbon abundance Xc(^^C) =0.20 for 
the M = 25Mq model. As shown in the table, mainly because of the difference in 
the ^^C(a,7)^^0 rate, the YU models presented in this paper have systematically 
smaller carbon abundance than the UN models. In Table II we show also final, Aff, 
He core, Mne, and CO core masses, Mqo- Here the latter two masses are defined 
where the mass fraction of H and He are less than 0.001, respectively. 
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Table I. Carbon mass fraction after core He 
burning for UN and YU models. 



UN model 


M{Mq) 


11.5 


12 


13 


15 


20 


22 


25 




0.36 


0.36 


0.33 


0.36 


0.34 


0.28 


0.26 


YU model 


M{Mq) 


11 


12 


13 


15 


18 


20 


25 




0.26 


0.25 


0.31 


0.24 


0.23 


0.23 


0.20 



Table II. Final, core and remnant masses for 
UN, YU and NH88 models. 



UN model 


M{Mq) 


11.5 


12 


13 


15 


20 


22 


25 


M{{Mq) 


11.2 


11.6 


12.7 


13.9 


17.5 


17.6 


17.7 


Mhc (Mq) 


2.8 


3.0 


3.3 


4.0 


5.3 


6.7 


7.9 


Mco (Mq) 


1.62 


1.72 


1.95 


2.51 


3.91 


5.42 


6.54 


Mrem (Mq) 


1.40 


1.43 


1.45 


1.51 


1.57 


1.66 


1.70 


Mg (Mo) 


1.26 


1.29 


1.30 


1.35 


1.40 


1.46 


1.50 


YU model 


M(Mq) 


10 


11 


12 


13 


15 


18 


20 


Mf (Mo) 


9.5 


10.5 


11.4 


12.1 


13.6 


16.2 


17.6 


Mho (Mo) 


2.6 


2.8 


3.1 


3.2 


4.2 


5.4 


6.2 


Mco (Mo) 


1.47 


1.61 


1.81 


1.87 


2.64 


3.58 


4.34 


M,em (Mq) 


1.29 


1.32 


1.44 


1.45 


1.64 


1.77 


1.96 


Mg (Mo) 


1.18 


1.19 


1.29 


1.30 


1.45 


1.55 


1.69 


NH88 model 


M(Mo) 


~ 13 


~ 15 


~ 20 


~ 25 








Mhc (Mo) 


3.3 


4 


6 


8 








M,em (Mq) 


1.27 


1.33 


1.61 


1.77 








Mg (Mo) 


1.15 


1.20 


1.43 


1.55 









3.1. Evolution of 10-13 Mq stars 

Here we describe the evolution of 10-13 Mq stars calculated by the YU code 
since these stars correspond roughly to the lightest Fe core collapse model, and the 
evolution after oxygen burning is quite different from more massive M > 13M0 stars. 
Table IIIII shows the differences in the way how the stars commence nuclear burning 
of Ne, O, and Si. In the table a letter 'C represents that the ignition occurs at 
the center of the core and 'off-C shows that the ignition occurs off-centrally. In the 
line of the composition, letters represent main nuclei at the center. The differences 
also appear in Fig. [3l in which evolutional tracks are shown in terms of the central 
temperature and density. In this mass range, a star forms a CO core around the 
Chandrasekhar limiting mass. The degeneracy of the core is the important property 
and affects its evolution. 
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Fig. 3. Evolution tracks of M = 9.82, 9.89, 10, 11, and 12Mq stars in terms of central temperature 
and density. 

We will present more detailed description in Takahashi, Yoshida & Umeda (2012, 
in preparation, see also Ref. [37|)). in which stellar evolutions are calculated in very 
fine grids. So far, we have found that M = 9.89M0 model ends up the Fe core 
formation, though M = 9.82Mq model does not form an Fe core. As described 
below the off-center Si-burning is more and more violent for less massive stars. In 
this mass range only 0.01 Mq difference in the zero-age main sequence mass would 
result in quite different evolution. However, the global evolutionary properties and 
final density structure of an Fe core is similar to the 10 Mq model presented here as 
far as an Fe core is formed. 

10 Mq star. For a 10 Mq star, the evolution from the zero age main sequence 
to the C burning phase goes roughly in the same way as more massive stars which 
end up with the Fe core formation. However, after an 0-Ne core formation through 
the C burning, different evolutionary aspects come about. 

Because the mass of the C-0 core of IAYMq is only slightly larger than the crit- 
ical mass for the Ne ignition, the 0-Ne core first contracts without nuclear burning. 
Eventually, the importance of the pressure of degenerate electron increases. Such 
semi-degenerate cores often show the inverse temperature distribution; temperature 
does not take its highest value at the center. This is due to the property of degen- 
erate electrons. Since neutrino cooling is more effective at higher density as well as 
at higher temperature, the 0-Ne core loses its heat by the cooling. As heat escapes, 
the core contracts and compressional heating supply energy into the gas of ions, 
electrons, and radiation. If the core is supported by pressure of ions, gas temper- 
ature increase as a result of compressional heating. On the other hand, if electron 
pressure supports the core, gas temperature is hard to increase because internal en- 
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Fig. 4. Mass fraction distribution for a lOM© star during a shell Ne+O burning, when logm pc 
7.9358. 



ergy of degenerate electrons does not depend on temperature but on density. What 
the high density 0-Ne core supports mainly is the pressure of degenerate electrons, 
so the denser the region is, the slower the temperature increase is. Owing to this 
temperature inversion, a shell Ne burning commences at an off-center region having 
higher temperature than the center. Since Ne and O have close ignition tempera- 
tures, once the shell Ne burning commences, heat generation increases temperature 
and the shell 0-ignition succeeds. Then a degenerate 0-Ne core is surrounded by a 
hot Si-shell (FigH]). These evolutionary properties are fully consistent with previous 
works (e.g., Refs. fT2|).l8])). 

Some of the nuclei of the Si cluster produced by the O burning can capture 
electrons in such a high density 0-Ne core. Since the core is mainly supported by 
the pressure of degenerate electrons, the decrease of the number of electrons causes 
the core contraction. Heat production by the contraction increases the temperature. 
When temperature at the base of the Si shell becomes high enough to ignite Si, the 
shell Si burning commences and transforms the Si shell into the Fe shell (Fig. [5|). 
The energy production rate of the burning from Ne to Fe is very high. So at every 
ignition, the core repeats the expansional cooling and compressional heating almost 
adiabaticallvl*^ and the burning front approaches the center little by little. After 
the burning front reaches the center of the core, an Fe core forms. Then the star 
collapse. Some vertical ascents in Fig. [3] show moments when a base of a off-center 
burning reaches the center of the core. For a IOMq star, the central temperature 
suddenly increase at log^Q pc = 8.4319 when the shell Ne-j-0-|-Si burning reaches the 
center. 



*•* This is seen in spikes around logj^g pc ~ 8 in Fig. |31 
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11 Mq star: For a 11 Mq star with a larger C-0 core of I.6IM0 than a WMq 
star, electron degeneracy of the core is less. This causes the difference in evolution 
after Ne and O ignite off-centrally. 

Because of the less degeneracy at the center, a difference between the maximum 
and the central temperature is less, so a shell 0+Ne ignition commences at deeper 
region than the WMq star (Fig. [6]). The temperature at the base of the shell 
increases by compressional heating as the same way as the IOMq star during the 
0+Ne burning front propagation. However, the burning front arrives at the center 
before the base temperature reaches the Si ignition and, thus, a Si core forms. This 
is the main difference between 10 and IIMq stars. 

Later, as the Si core contracts, outer 0+Ne shell burning repeats moving its 
base outward, while neutron rich nuclei such as ^^S and ^''Ti are produced in inner 
region. Degeneracy of electrons also increases in the contracting core and the above- 
mentioned temperature inversion appears. When the central density reaches the 
value of log]^o PC = 8.9867, the first off-center Si ignition occurs at Mr = O.5933M0. 
This ignition causes the core to expand adiabatically. After the burning terminates, 
the core contracts adiabatically again. When the central density becomes log^g pc = 
9.0253, the second off-center Si ignition occurs at Mr = O.11O3M0 (Fig. [7]). After 
the second off-center Si ignition, the burning front reaches the center. The core 
transforms into an Fe core, then collapses. 

12 Mq star. A I2M0 star evolves similar to a UMq star but the properties 
of Ne and O ignitions are different. Because it has a larger C-0 core of 1.85Mq, 
temperature takes its highest value at the center of the 0-Ne core. Therefore, Ne 
and O ignite at center. Note that Table Hill shows that when O ignites at the center. 
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Fig. 6. Mass fraction distribution for a IIM© star during a shell Ne+O burning, when logjQ pc 
7.6181. 
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Fig. 7. Mass fraction distribution for a IIM0 star during a shell Si burning, when logjQ pc ~ 8.0481. 



central composition is 0-Si, and not 0-Ne. This is just because Ne has burned ahead 
of O, producing a 0-Si core. 

After Si core formation, the evolution proceeds as the same way as a UMq star. 
Neutron rich nuclei are produced around the center after several shell Ne+O burning 
at outer regions of the contracting core. Then off-center Si burning occurs (Fig. [5]), 
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MJMq 

Fig. 8. Mass fraction distribution for a 12Mq star during a shell Si burning, when log^Q pc = 7.9694. 
and after the front reaches the center an Fe core forms. 



Table III. Ignition properties of 10 — ISM© 
stars. 



C ignition 


lOMo 


IIM0 


12A/0 


13M0 


at the center or not 


C 


C 


C 


C 


central composition 


C-0 


C-0 


C-O 


C-O 


Ne ignition 


1OM0 


11M0 


12Mq 


13Af0 


at the center or not 


off-C 


off-C 


c 


C 


central composition 


0-Ne 


O-Ne 


O-Ne 


O-Ne 


ignition 


lOAf0 


IIMq 


I2M0 


13Mq 


at the center or not 


off-C 


off-C 


C 


C 


central composition 


0-Ne 


O-Ne 


0-Si 


O-Si 


Si ignition 


lOAfo 


llA'/0 


12M0 


13M0 


at the center or not 


off-C 


off-C 


off-C 


off-C 


central composition 


0-Ne 


Fe-Si 


Fe-Si 


Fe-Si 



3.2. Final mass and metallicity 

In Fig. [9] we show the metalhcity dependence of the final and core masses of 
the YU models.!^ We see clear metallicity dependence in the final mass among the 
stars with M > 20Mq. In Z = 0.02, M = SOM© star has the maximum final mass: 
Mf = 21.3M0. M > 4OM0 stars indicate a roughly constant final mass of ~ WMq. 
These stars lose all of H and He layers and become Wolf-Rayet stars. Effect of 
mass loss in Z = 0.01 stars is less than that in Z = 0.02 stars. Increasing the main- 
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Fig. 9. The final mass (panel (a)), the He core mass (panel (6)), and the CO core mass 
(panel(c)) with the relation to the main-sequence mass. The adopted metallicities are Z — 
10"*, 0.001, 0.004,0.01, and 0.02. In panel (&), stars more massive than the attached number of 
the main-sequence mass (Mq) evolve to WO stars. The surface He mass fraction of the stars is 
about 0.2. 
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sequence mass, the final mass also increases in M < 30Mq and 40 < M < 80Mq but 
it decreases in 30 < M < 4OM0 and M > 80Mq. Metal poorer stars have similar M 
dependence of the final mass. In this figure the zero-age main-sequence mass where 
the mass loss becomes effective and the final mass for a given main-sequence mass 
become larger for metal poorer stars. We do not see clear mass loss effect among 
M < 80Mq stars with Z = W'^. 

The metallicity dependence of the He core is small (Fig. El^fc)). On the other 
hand, most of the He layer has lost during the evolution in M > 40, 50, and 100 
Mq stars for Z = 0.02,0.01, and 0.004. These stars lose the He/C/0 envelope and 
their surface He abundance decreases during the He burning. They become C- and 
0-enriched Wolf-Rayet stars (WO stars) with the surface He mass fraction of ~ 0.2. 
The mass of the He-rich shell is 1.4 — 5.5Mq in Z = 10~^ stars, in which mass loss 
effect is small. Since the fraction of the He layer is small, mass loss brings about the 
removal of He layer rather than the reduction of the He core mass. 

The metallicity dependence of the CO core mass is seen in Z > 0.004 and 
M > 40 — 50Mq in Fig. [O^c). The CO core mass is roughly constant with ~ WMq 
in Z = 0.02 and M > 4OM0 stars. In Z = 0.01, the maximum mass of the CO core 
is 22Mq of M = 80Mq star. On the other hand, the CO core mass can be more 
massive in metal-poorer stars. The CO core mass is 35 — 40Mq in Z < 0.004 and 
M = IOOMq. 

§4. Nucleosynthesis 

Nucleosynthesis in these massive stars occurs mainly in two stages. One is 
before core-collapse and the other is during supernova explosion. It is well known 
that inside these massive stars nuclear fusion takes place up to Fe synthesis. Before 
gravitational core collapse, Fe core is produced at center and lighter elements form 
onion like structures from inside to the surface (Fig. [TO]l . 

Nucleosynthesis during the supernova explosion is usually called explosive nucle- 
osynthesis. During the explosion, shock wave propagates out of Fe-core to the stellar 
surface. Behind the shock wave matter heats up and nuclear burning takes place. 
Explosive nucleosynthesis is important for the synthesis of Si and heavier elements 
(e.g., Ref.lE])). 

In our previous works (e.g., Refs. [T6j) . [22j) . rr7|) . l38]) ) we use a simple ID model 
for supernova explosion to calculate explosive nucleosynthesis. We inject thermal or 
kinetic energy below the mass cut or just above the Fe core to initiate supernova 
shock. We may call this procedure as 'instant energy injection'. 

We use a ID PPM code for the hydrodynamical calculations and solve small 
alpha- network together to calculate nuclear energy generation. Then we calculate 
detailed nucleosynthesis post-processingly by solving a large nuclear reaction net- 
work. 

Explosive nucleosynthesis should depend on how the star explodes, but we do 
not know yet how gravitational collapse leads to the explosion. Therefore, there are 
still several proposals for the successful supernova explosions. Fortunately the prop- 
erties of supernova shock outside the Fe core is almost independent of the explosion 
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Fig. 10. Mass fraction distribution of the supernova progenitors evolved from 15 (a) and 20 Mq 
(6) stars with the metallicity of Z = 0.02. The corresponding final masses are 13.6 and 17.6 
Mq. The mass fraction distribution in the H-rich envelope outside the range of each figure is 
the same as that at the outer edge of the figure. 



mechanism, because in most models the supernova shock gains energy inside the Fe 
core. Then, as far as spherical symmetry is assumed, explosion energy E is the only 
parameter which determines the properties of the supernova shock. For example, 
kinetic to thermal energy ratio does not much affect the propagation of the shock 
outside the Fe core, because it quickly converges into the same solution. In this case 
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'instant energy injection' method is sufficient to calculate nucleosynthesis. For most 
cases, nucleosynthesis of Fe peak elements including Zn and lighter elements can be 
safely calculated with this method. In Section 7, We will discuss the limitations of 
this simple procedure. 

Behind the shock front, temperature distribution is almost homogeneous, and 
the radiation energy is strongly dominated. Then the radiation temperature T can 
be related to the explosion energy, E,hY E = (47r/3)ar3r^, where a is the radiation 
constant and r is the shock radius. Thus the maximum temperature after shock pas- 
sage, Ts, for the mass elements located at radius R in the progenitor is approximately 



For a given progenitor model and explosion energy E, the propagation of shock 
and thus the time evolution of density distribution is also determined uniquely. Be- 
hind the shock front, pressure and entropy are also radiation dominated and can 
be written as = aT'^/3 and = (4/3)aT^/p. Since the shocked region cools 
roughly adiabatically at first, the radiation entropy is often convenient to specify the 
explosion. 

In the inner most region of supernova ejecta where Tg is higher than 5x 10^ (K), Si 
burns completely. In such a region, at first Si mostly decomposes to alpha-particles, 
and this is endothermic reaction. Then as the temperature decreases alpha particle 
starts to recombine and alpha-rich freezeout nucleosynthesis takes place. This phase 
is exothermic. The energy changes due to these processes are typically about ten 
percent, and thus above equation ()4-ip still provides a useful approximation. 

In this 'complete Si burning' region ^^Ni is dominantly produced but also Co, Ni, 
and Zn are mostly produced here. The final abundance depends on entropy because 
the mass fraction of alpha-particles during the alpha-rich freezeout phase increases 
with entropy.--' 

For example, one of the authors has shown that high entropy during the hyper- 
nova explosion explains well the larger ratio of [Zn/Fe] and_[Co/Fe] towards lower 
metalhcity, [Fe/H], in extremely metal poor (EMP) stars-H^'El®,!!! They also 
discussed that simple ID model cannot explain the abundance of EMP stars by the 
hypernova models, because Fe or ^^Ni is over-ejected to satisfy the large [Zn/Fe] 
ratio. This problem and the solution is described in the subsection 4.3. 

When Ts is between 4 — 5 x 10^ (K), Si burns incompletely. This region is called 
incomplete Si burning region. Main products here is Cr, Mn, and ^^Ni. For Tg ~ 
3 — 4 X 10^ (K), oxygen burns partially and Si, S, Ar, Ca, Ti, V are main products. 
In these regions also nuclear burning produces some amount of energy. In summary, 
previous works have shown that explosive nucleosynthesis up to Zn can be described 
well by instant energy injection models, though hypernova models over-produce Fe 
in simple spherically symmetric models. 

So far we have not mentioned the effects of neutrino process (I'-process), but the 



process is important for some elements such as Li, B, F, and Mn.SD -US. 33 'HI -113.133,111 ,113 ,111 ,111, El 



Of course the neutrino emission depends on the explosion model and thus we have to 
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assume a model to include the effects (see Ref . [38|1 for the detail) . Fig. [12] shows the 
mass fraction distributions of Li, B, F, Sc, V, and Mn of the supernova ejecta evolved 
from a 15 Mq star with Z = 0.02. In this calculation, we assume that the neutrino 
luminosity decreases exponentially with the decay time of 3 s. The total neutrino 
energy is set to be 3 x 10^^ erg. The neutrino spectra are assumed to obey Fermi 
distribution with zero cliemica,! pot6iiti3<l. TI16 temperatures of ^/j,,t; ^/j.,t 

and Ue, Ve 
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MJMq 

Fig. 12. Mass fraction distribution of Li, B, F, Sc, V, and Mn produced through the i^-process 
in the supernova ejecta evolved from a 15 Mq star with the metallicity of Z — 0.02. Solid 
and dotted hnes indicate the mass fractions of the elements with and without the i^-process, 
respectively. 

are set to be 6 MeV/k and 4 MeV/k, respectively, as in Ref. [38]) . Almost all Li and 
B are produced through the z/-process in the 0/Ne-layer and He-layer. When the 
z^-process is not taken into account, the mass fractions of Li and B are smaller than 
the lowest value in this figure. Large F production through the i^-process is obtained 
in the 0/Ne layer (see also Refs. H^ .ISU j) ). A part of F is produced in the explosive 
He burning from ^^N. Although the F production from ^^N strongly depends on the 
metallicity, the F production through the z^-process scarcely depends on the metal- 
licity and, thus, it is important throughout the Galactic chemical evolution. Parts 
of Sc, V, and Mn are produced though the i^-process, especially in the complete and 
incomplete Si-burning regions. 

4.1. Xc(^'^C) and abundance pattern of Ne to Ca 

In the previous section, we mentioned that the carbon abundance after helium 
burning is important for the abundance of Ne to Ca. Here we look the results of 
UN and YU models more closely. In Fig. [T3]we show [Xj/0] vs. atomic number. 
Here Xi represents Salpeter's initial mass function (IMF) weighted yields of Ne to 
Ca isotopes, and [A/B]= logio(YA/^B) — logio(^A/yB)05 where Ya and 1b are the 
abundances of elements A and B. If any point is far from the solar value, [Xj/O]=0, 
the model fails to explain the present day abundance, though chemical evolution 
model should be applied for a detailed discussion. This figure show that for both 
models, each point is roughly in the it 0.3 range, thus these models would yield 
roughly solar abundance pattern. 

Since UN model have larger Xq{^'^C) after helium burning, this model has larger 
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Fig. 13. The IMF weighted yield ratios [Xi/0] as a function of atomic weight. Red and blue lines 
indicate the yield ratios in YU model and UN model, respectively. 



Ne/0 and Na/0 ratios than the YU model as expected. Larger Xc(^'^C) means that 
smaller Xc(^^0) and thus explosive oxygen burning products, such as S, Ar, Ca, 
would be less abundant. This is also seen in the figure. 

The UN yields calculated with the same parameter choices was applied to the 
Galactic chemical evolution model in Ref. [5T]) . It was shown that the model could 
reproduce the observed abundance pattern reasonably well. If we look the results 
closely, however, [Na/Fe] is slightly larger than zero at [Fe/H]=0, and [Ca/Fe] and 
[Ar/Fe] are slightly smaller than zero at [Fe/H]=0. This suggests that slightly larger 
^^C(a,7)^^0 rate may reproduce the solar abundance better at [Fe/H]=0. Using 
the YU codes we will investigate the best choice of the ^^C(a,7)^^0 rate, or the 
Xc(^^C), in a forthcoming paper. 

4.2. Zn abundance in EMP stars and asphericity of hypernovae 

Observations of extremely metal poor (EMP) stars, which are defined as halo 
stars with [Fe/H] < —4, revealed that [Zn/Fe] and [Co/Fe] are larger than higher 
metallicity stars (e.g., Ref. [52]) ). Ref. \W\i proposed that these large ratios are ex- 
plained by the ejecta of complete Si burning of hypernova. This idea explains also 
why such EMP stars have low [Fe/H]. In the early universe where the metal abun- 
dance in inter stellar matter (ISM) was low, the Fe to H ratio was determined by the 
amount of Fe produced by a supernova over that of hydrogen swept by the super- 
nova shock. Since the swept hydrogen mass is roughly proportional to the explosion 
energy, the ejecta of a hypernova tends to have lower [Fe/H] than that of a supernova. 

As mentioned above, however, there is one problem for this idea. Zn and Co are 
produced mainly in the complete Si burning while Fe or ^^Ni is produced both in 
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complete and incomplete Si burning regions. Therefore, [Zn, Co/Fe] becomes larger 
when relatively more material is ejected from the complete Si-burning region. Under 
the assumption of spherical symmetric explosion, this means that the mass-cut has 
to be taken sufficiently deep, then the ejected mass of Fe or ^^Ni also increases. 
On the other hand if too much Fe is ejected 0/Fe ratio, for example, becomes too 
small compared with observations. To quantify this problem Ref. I16p introduced the 
'mixing- fallback model'. In the model, yields are calculated as follow. First, inner 
most matter is mixed between the mass coordinate Mjn and Mout- For hypernova 
models Mout is set to the upper boundary of the incomplete Si burning region, and 
Min is chosen sufficiently deep to eject Zn. Then it is assumed that only the fraction 
/ of the matter is ejected, or the fraction 1 — / fallbacks, from this region. To explain 
the abundance pattern of EMP stars / ~ 0.1 is required. 

For supernova models with M < 25Mq, f=l is fine to reproduce abundance 
pattern of very metal poor stars with —3 < [Fe/H] < —2.'^ This suggests that the 
explosion mechanism of supernova and hypernova is different. 

Interestingly almost the same yield as the / ~ 0.1 hypernova yield could be 
obtained when we calculate jet-like hypernova models in 2D.^^' In such models, it 
is assumed that an unknown central engine ejects jets along the polar directions. 
These jets can blow up the entire star above the Si layer, but the explosive Si 
burning mainly takes place along the jet directions. As a result the mass fraction 
of complete Si burning products becomes smaller than the spherically symmetric 
models. In this paper also in §6.3 we show such a 2D calculations applied to an 
initially M = IIOMq star. Since observed hypernovae show some asphericity,'^ 
this jet-like explosion model is certainly interesting, though we need to know the 
mechanism for the central engine. 

It is sometimes claimed that Zn and Co in EMP stars may be explained by the 
innermost matter from 'hot bubble' region of normal supernovae (e.g., Ref. "^5)). 
Ref. i56l) explored this possibility and discussed that this explanation is not likely 
because fine tuning of Ye{= 0.500 — 0.501) is required for a 0.06 Mq ejecta to produce 
sufficient amount of Zn and Co. These conditions are difficult to satisfy for the hot 
bubble matter of supernovae, while it is relatively easy for the innermost ejecta of 
hypernovae. 

4.3. Nucleosynthesis of (weak) r-process elements 

It is quite likely that during explosion supernovae produce elements heavier 
than Zn, because some EMP stars have r- and weak-r process elements. For exam- 
ple CS22892-052 shows r-process abundance pattern which is roughly same as the 
present-day r-process pattern.'SS This suggests that the r-process nucleosynthesis 
is almost universal. There are other classes of stars which show more abundance 
of weak r-process elements such as Sr, Y, Zr than the universal r-process pattern. 
There are also some stars having enhancement of Mo and Ru as well as Sr, Y, Zr 
stars 

As mentioned above spherically symmetric instant energy injection models can- 
not synthesize such elements because to synthesize such elements much larger entropy 
and/or more neutron rich environment are required. On the other hand Fe-peak ele- 
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ments including Zn can be synthesized from the matter with Ye ~ 0.50. Most recent 
ID core collapse supernovae simulations have shown that the innermost matters in 
the ejecta have Ye — 0.50 or even greater than 0.5 because of the interaction with 
neutrinos (e.g., Refs. [MDEQDISI])). 

Ref. HOj) and l56p explored the conditions for synthesizing weak r-process elements, 
Sr, Y, and Zr. Typically innermost layers of normal supernova explosions have 
entropy per baryon s/A; ~ 5 and those of hypernovae have s/k ~ 15. With such 
entropy weak r-process elements are not produced if 1^ ~ 0.50. Thus they relaxed 
the constraint Ye — 0.50 and considered Ye as low as 0.45. Such low Ye matter cannot 
be ejected as long as we consider exact spherically symmetric explosion. However, as 
shown in Ref l62|) in 2D calculations inner materials are mixed by convection during 
explosion and small amount of low Ye ~ 0.45 and higher entropy s/k ~ 40 — 50 
matter, may be ejected. Ref. |40j) showed that s/k > 15 model could produce Sr, Y, 
Zr explosively. Therefore, such low Ye matter ejected due to the multi-dimensional 
effects can explain the origin of weak r-process elements. 

We note that s/k < 50 models could produce up to Zr but not heavier elements. 
Because of the lack of observations it is not clear yet weak r- stars which have 
enhancement of Sr, Y, Zr always have Nb-Mo enhancements as well. If this is the 
case, innermost matter of hypernovae or hot bubble of normal supernovae considered 
in Ref. H0|) may not be the main site for the weak r-process synthesis. To produce 
Nb-Mo much larger entropy {s/k ~ 150) is required.!^ It is currently not clear 
which astronomical site has such an environment. 

§5. Remnant neutron star mass 

In this section we present remnant neutron star mass in both the UN and YU 
models and discuss their implications. There are several factors to determine the 
mass such as the CO-core mass, explosion energy and EOS. Among them the most 
important factor is the CO-core mass, which is determined by the stellar evolution 
before core-collapse. Larger CO-core leads larger remnant mass because it typically 
leads larger Fe-core and, more importantly it increases the amount of mass above 
the Fe-core. 

Once pre-explosion density structure is given, explosion energy determines the 
mass-cut or remnant neutron star mass Mj-em- For a given progenitor model, larger 
explosion energy blows up more materials above the Fe-core leading smaller mass- 
cut. In this paper, however, as well as our previous works, we do not determine 
mass-cut in this dynamical way but determine it by the amount of the ejected ^^Ni 
amount, M(^^Ni). This is because M(^^Ni) is rather sensitive to explosion energy 
when we determine the mass-cut dynamically. In reality each SN may eject somewhat 
different amount of ^^Ni. However, when we apply the supernova yields for example 
to the Galactic chemical evolution it is more useful to provide averaged yields. For 
this purpose it is better to determine the mass-cut by M(^^Ni). 

We set the explosion energy as E^xp = 10^^ erg, and the ejected ^^Ni amount as 
0.07 Mq to determine M,em- M{^^m)=0.07 Mq is the value for the SNIOSTA'^^' and 
considered typical value for normal core-collapse SNe. Strictly speaking to eject the 
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same amount of Ni more massive core requires larger E^xp- However, the resultant 
Mi-cm is not much different, so we fix the -Egxp for simplicity. 

In Table II, we show the baryon mass of the remnant neutron star mass by Mi-em- 
It correlates with the CO-core mass, which in turn correlates with the He-core mass. 
The UN model leaves smaller remnant than the YU model for the same initial mass, 
M, mainly because Xc(^^C) after He-burning is larger and thus Mco is smaller as 
described in §3. 

As described, the choice M{^^Wi)=0.07 Mq is a typical value, but this is of course 
not the unique vale (see e.g., a review in Ref. lEj)) Therefore, we should remind this 
fact when we compare with observations. For example, when M(^^Ni)=0.01 Mq in 
the YU model, M^em = (1-38, 1.41, 1.53, 1.53, .172, 1.86, 2.04) Mq for M = (10, 11, 
12, 13, 15, 18, 20) Mq, respectively. 

The Mrem is the baryon mass and is not the observable neutron star mass. 
The observed mass is about 10 percent smaller due to the general relativistic effect 
and called gravitational mass (see e.g., Ref. I65|) ). The conversion from baryon to 
gravitational mass depends on the EOS of nuclear matter. In Table II we show the 
gravitational mass, Mg, for the UU model in Ref. i66p using UV14+UVII EOS^ for 
nuclear matter. 

5.1. Comparison with the observed mass 

Here we compare the Mg of our models with observed neutron star (NS) masses. 
NS masses are most accurately determined when they are in double NS systems. 
In this case if two or more post-Keplerian parameters are obtained, NS masses are 
precisely determined (e.g., Ref. [68 ]) .l69 p ). In Table HI we show the NS masses and 
errors for these stars. Interestingly the mass distribution has peak at 1.33 Mq with a 
small dispersion 0.06 Mq. These stars may be considered to keep their birth mass,'^ 
so they are suitable to compare with our results. NSs in binaries with high mass 
companions are also considered to roughly keep their birth mass. ReL^^E^ discussed 
that the most likely values of the central mass and dispersions for these NSs are 1.28 
Mq and 0.24 Mq, respectively. 

First we note that lowest NS mass 
in the table is Mg = 1.16 - I.2OM0, 
and this is consistent with the lowest Mg 
model in the YU model M = lOM©. As 
described above this model may roughly 
correspond to the lightest Fe-core col- 
lapse supernovae. For a smaller mass 
star such as M = 9.5 Mq or less, the 
star once forms a 0-Ne degenerate core 
and its collapse may lead a weak super- 
nova explosion .'SI 

In the calculation by Ref. I14|) . the 
remnant neutron star has Mrem = 
1.36M0 or Mg ~ 1.23M0. This seems 
to be larger than the observed minimum 



Table IV. NS masses in double NS systems 
and errors. See references for Ref. I70II and 
therein. 



Name 


Mass {Mq) 


Error (Af( 


J0737-3039 


1.3381 


0.0007 


pulsar B 


1.2489 


0.0007 


B1534+12 


1.3332 


0.0010 


companion 


1.3452 


0.0010 


J1756-2251 


1.40 


0.02 


companion 


1.18 


0.02 


J1906+0746 


1.248 


0.018 


companion 


1.365 


0.018 


B1913+16 


1.4398 


0.002 


companion 


1.3886 


0.002 


B2127+11C 


1.358 


0.010 


companion 


1.354 


0.010 
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mass NS, though these 0-Ne supernovae 
have not been studied well and currently 

only one progenitor model exists pD''^ therefore we do not know the general prop- 
erties of these SNe yet. 

For the observed double NS typical mass range is Mg = 1.27 — I.SOMq as 
mentioned above. This corresponds to M ~ 11.5 — 20Mq for UN model and M ~ 
12 - UMq for YU model. 

Largest NS mass in the table is Mg = 1.44Mq and this corresponds to M ~ 
22M0 for the UN model, and M ~ 15Mq for the YU model. It is usually said that 
the border of NS and blackhole formation is at around M = 20 — 25Mq. Therefore 
it appears that the Mg range of the UN model is consistent with the double NS mass 
range, though we have not confirmed yet if the UN model with M < UMq would 
result in Mg ~ I.2M0 NSs. 

If we adopt the YU model, the mass range of double NS corresponds to rather 
narrow range, M ~ 12 — 14M0, and it is difficult to understand the reason why. 
However at this moment we cannot say that the UN model represents the reality 
better, because the number of the double NS systems are still limited. In Ref. [70|1 
using the model of Ref. I73p they also discussed that the narrow mass range of Mg in 
double NS systems are difficult to understand and suggested a particular and rare 
formation channel for the systems. On the other hand, the NS mass in binaries with 
high mass companion ranges Mg = 1.04 — 1.52M0.'2SI Although this estimate is more 
uncertain than that in double NS systems, this wider range would be more easily 
understood in the YU model. 

In summary NS mass observations would certainly provide interesting and im- 
portant information to constrain the progenitor model, Xc(^^C) and ^^C(a,7)^^0 
rate. Interestingly all the observed NSs mentioned above have relatively small mass 
compared with the "recycled" NSs. NSs with white dwarf companions, millisecond 
pulsars and in low-mass X-ray binaries are called recycled and currently undergoing 
accretion. It is now known that these NS can surely be as massive as Mg = 2.OM0 as 
the case for J1614-2230.'^ Observations and theory are consistent in the sense that 
such massive NSs are rarely formed by the normal supernova explosions. Theoreti- 
cally, when we fix the explosion energy, iVfj-gj^ increases rapidly around ikf ^ 25A^0. 
This explains why massive NSs are rarely formed. It will be very interesting to 
constrain observationally that how massive NS can be formed at birth. 

§6. Other recent works of our group 

In this paper we have described our recent work on stellar evolution with some 
new results. Here we briefly review other works of our group not mentioned above. 

6.1. (Magneto- )Hydrodynamical Simulations 

T. Kuroda and H.U. have developed a 3D magnet-hydrodynamical general rel- 
ativistic code with adaptive mesh refinements.'^ This code was used to follow 
gravitational Fe core collapse and to calculate spectra of gravitational waves. We 
will apply this code to various progenitor models to explore explosion mechanism 
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and nucleosynthesis. 

6.2. Special Relativistic Hydrodynamical Simulations for GRB Jets 

S. Okita and H.U. have developed a 2D special relativistic hydrodynamical code. 
Using this code, Okita and H.U.'^SI' explored the conditions for successful ejection of 
ultra-relativistic jets in the collapsar model of GRBs. This code is also applied to 
explore various aspects of explosions and nucleosynthesis in SNe and GRBs. One 
example is given in the next subsection. 

6.3. Evolution and nucleosynthesis in very massive stars which end up Fe core- 
collapse SNe 

Motivated by the discoveries of very luminous SNe, such as SN 1999as (SN Icf^ 
and SN 2006gy (SN IIn),E3 H.U. and K. Nomoto^^' calculated evolution of metal- 
poor massive stars in the mass range M = 20- IOOM0 with metahicity Z = 10"^ to 
study how much ^^Ni is produced in core collapse SNe (CCSNe). They found that 
^^Ni of ~ 13Mq can be produced for the 100 Mq star. This amount is sufficiently 
large to explain SN 1999as if the SN shines by the ^^Ni decay. 

Since the actual bright SNe Ic usually appears in a host galaxy having metal- 
licity larger than Z = 10~'^, mass-loss effect becomes important. Therefore, T. Y. 
& H.U.'^S studied the uncertainties in mass-loss in detail for Z = 0.004 which cor- 
responds to the host galaxy of SN 2007bi. SN 2007bi was a very bright SN Ic and 
3.6 — 7.4M0 of ^^Ni is required if it shines by radioactive ^^Ni. First Gal- Yam et 
al.ESJ discussed that the SN was a pair-instability SN (PISN) but Moriya et al.^J 
showed that they can be a CCSN if a ~ 43M0 CO star explodes with £^exp = 3 x 10^^ 
erg. 

Since PISN model is much more massive than the CCSN model, these mod- 
els predicts quite different light curves. Unfortunately without the LC data well 
before the maximum light one could not distinguish these two models. T.Y. and 
H.U.ISSJ found that for Z = 0.004, required ZAMS ranges are M = 110 - 28OM0 and 
M = 515 — 575M0 for CCSN and PISN model, respectively. They provided that, if 
the progenitor was a single star and assuming the Salpeter's IMF, the ratio of the 
probabilities of CCSN to that of PISN appropriate for SN 2007bi is 42. We should 
remind, though, for the CCSN model we are assuming that such a large CO star can 
explode energetically. So far no one could have shown such a explosion from the first 
principle calculations. 

Under the assumption that it can explode Okita, H. U., and T. Y. simulated 
spherically symmetric and axisymmetric jet-like core-collapse supernova explosions 
of a Mf = 43.2M0 WO star in Ref. |20D.ED'[S1 This progenitor is evolved from a 
M = IWMq star with the metalhcity oi Z = 0.004. The CO core mass is Mqo = 
38.2M0. For spherically symmetric calculations the same method as described in §4 
is used. For the jet-like explosion the same code mentioned in §6.2 was used. After 
the explosion simulations, nucleosynthesis is calculated post-processingly. First we 
investigate the explosion-energy dependence of the ^^Ni amount. We found that ^^Ni 
larger than 3 Mq, enough to reproduce the light curve of SN 2007bi, is produced 
in the supernova ejecta if the explosion energy is larger than 2 x 10^^ erg. The 
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Fig. 14. The relation between the ejected ^^Ni amount and opening angle ^op in axisymmetric 
supernova model. The progenitor is a 43.2 Mq WO star. The explosion energy is set to be 
3 X 10^^ erg. 



investigation by the jet-like explosions with different opening angles ^op with i?exp = 
3 X 10^^ erg indicated that the ejected ^^Ni yield strongly depends on the opening 
angle. Fig. [T3] shows the relation between the amount of the ejected ^^Ni and the 
opening angle. Although the spherical explosion releases the amount of 4 Mq ^^Ni 
in the ejecta, the ejected ^^Ni is much smaller when the opening angle is smaller 
than 78°. The ejected amount of ^^Ni is smaller than 2.3 Mq for the opening angle 
smaller than 68°. This suggest that if SN 2007bi was a CCSN, it exploded with large 
opening angle or the explosion was more energetic than 3 x 10^^ erg. 

6.4. Pop III very massive stars: Evolution, Dark stars, Stability 

The formation of first generation or Population (Pop) III stars in the Universe is 
considered to be quite different from that of later generation stars. This is because 
when the first generation stars were formed there was no metal yet, thus radiative 
feedback from the proto star was weak and might not be able to stop gas accretion. 
Also gas accretion rate itself was larger than the present universe. As a result it was 
expected that typical mass of Pop III stars were 100 Mq or more (e.g., Ref. I83|) ). 

In Ohkubo, H.U. et al.j'SS they calculated evolution and nucleosynthesis in the 
very massive 500 and 1000 Mq stars assuming that they explode energetically. Next, 
by taking account of a realistic mass accretion rate, Ref. [85j) calculated Pop III star 
evolution with the mass accretion. They found that typically 500 — IOOOMq stars 
were formed. It is notable that this mass range is larger than the often said mass 
range for a Pop III star, M = 140 — 3OOM0. If the mass range is in the 140 — 3OOM0 
range, most first stars should have exploded as PISNe. However, the abundance 
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patterns of EMP and HMP stars with [Fe/H] < —4 are not consistent with the 
yields of PISNe.liS 

These works, however, neglected the effects on the accretion disk likely formed 
around the first stars. Ref. [86]) and 1871) investigated such effects and showed that 
the accretion disk was evaporated by the radiative feedback. Then the accretion 
onto a first star may stop before its mass grows more than 100 Mq, preventing the 
formation of a PISN.E3 

Even though it is true, there is still a possibility for the formation of 1000 Mq 
stars. First stars are formed around the center of a dark halo where the density of 
dark matter is much higher than other places. Refs. [88 |) .l89 p showed that if the dark 
matter is self-annihilating WIMPs the self-annihilation energy is sufficiently large to 
sustain the first stars. Such stars are called "dark stars". These stars typically have 
much larger radius, thus surface temperatures are lower and radiative feedback is 
expected to be weak. Hirano, H.U. and N. YoshidsP^ calculated the dark star model 
with mass accretion. We confirmed that the stars keep large radius until it reaches 
the main sequence stage at around M ~ lOOOM©. 

H.U. et alW^ also considered the accreting dark star model including the effects 
of captured dark matter annihilation that becomes important after the main sequence 
stage. They showed that if the baryon-dark matter scattering cross section is as 
large as cr = 10~^^cm~^, the dark stars could grow more: it could be lO^tolO^ 
MQdepending on the mass accretion rates. 

It has been known that such huge stars are vibrationally unstable against epsilon 
mechanism (e.g., Ref. 92)). As far as the star is sustained by the DM annihilation 
and not nuclear burning, the epsilon mechanism does not operate. However after 
the main sequence stage, nuclear burning soon dominates unless the captured dark 
matter effect is important. Sonoi and H.UJ93P explored the stability of such very 
massive stars against the epsilon mechanism. They found that the amount of mass 
loss is less than 10 percent of the whole stellar mass. 

6.5. Dust Formation in Supernovae 

H.U. has been working also on dust formation in supernovae with T. Nozawa, T. 
Kozasa and collaborators. For example, Nozawa et al.'® calculated the dust yield 
for Pop III supernovae including CCSNe and PISNe, and showed that large amount 
of dust grains would be produced in the early universe by these SNe. Ref. 
calculated the dust destruction by the reverse shock in a supernova remnant. The 
theory of dust formation was also applied to actual supernovae. Ref. [96]) and [97|) 
for SN 2006jc and Cas A supernova remnant, respectively, showed that the observed 
data are reasonably reproduced by the theory. It was also applied to SNe laW^ These 
dust grains in supernovae play important roles in the star and galaxy formation, and 
chemical evolution in the galaxies. 

6.6. Presolar Grains from Supernovae 

Presolar grains are recovered from primitive meteorites or interplanetary dusts 
and are identified as the grains having very large isotopic anomalies compared with 
the solar-system materials (e.g., reviews [99j) . ll00p ). Observed isotopic ratios of the 
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grains are considered to indicate the traces of the nucleosynthesis in stars at their 
birth or Galactic chemical evolution. Small amount of presolar grains are considered 
to originate from supernovae. They mainly indicate the excesses of ^^C, ^^N, and 
^^Si. Some supernova grains show evidence for original presence of radioactive *^^Ti 
in Ca isotopic ratios, which strongly supports their origin. Isotopic ratios of heavy 
elements such as Mo and Ba have been also observed. However, it is still difficult 
to reproduce observed isotopic ratios by supernova models. Bulk composition of 
supernova ejecta of supernova models does not reproduce observed isotopic ratios of 
supernova grains. ^^^^ In order to reproduce observed isotopic ratios, inhomogeneous 
mixing is required. 

TY and HashimotcP^ and TYl^^ investigated supernova mixtures reproducing 
several isotopic ratios of supernova originating SiC and graphite grains. They di- 
vided supernova ejecta into seven different layers and investigated the mixing ratios 
reproducing C, N, O, Al, Si, and Ti isotopic ratios as many as possible for individ- 
ual grains. The mixing ratios of the mixtures strongly depend on the reproduced 
isotopic ratios. The main component of the mixtures are the innermost Ni layer 
and the outer He/C and He/N layers, so that inhomogeneous mixing of supernova 
eject is necessary. TY, HU, and NomotcP® investigated the supernova mixtures of 
different stellar masses reproducing Si isotopic ratios of supernova grains. They ob- 
tained that less massive supernovae with M < 15Mq and hypernovae are preferable 
to reproducing Si isotopic ratios of supernova grains. 

§7. Discussions and Future Prospects 

7.1. Massive star evolution 

In this paper we explained the differences in the newly developed efficient YU 
code and previous UN code, and shown that the YU code yields reasonable results 
as shown in §3 in some detail. We need such an efficient code because the study of 
massive star evolution still requires heavy amount of calculations as described in the 
followings. 

As shown in Tables I and II, the results of UN and YU models presented in 
this paper are different mainly because Xc(^^C) is different. The values of Xc(^^C) 
sensitively depend on the ^■^C(a,7)^^0 rate and treatment of convection. Since 
it is currently not possible to determine the ^^C(q;,7)^^0 rate, we need calculate 
for various choices of the rate to find a best set to fit observations. Traditionally 
nucleosynthetic argument shown in §4.1 has given the most stringent constraint on 
Xc(^^C) and the ^^C(a,7)^^0 rate. However, as shown in Fig. 11, it is not easy to 
distinguish even the UN and YU models because abundance data are only given for 
the IMF weighted yields. 

We propose here that NS mass distribution will give alternative and independent 
constraints on the Xc(^^C) and -^^C(a,7)^^0 rate. As shown in §5 and Table II, UN 
and YU models predict different mass distribution for the remnant NSs. Though 
the differences are not so large, observed NS masses are given in very fine precision 
for the double NS systems (Table IV). Therefore, by calculating stellar evolution in 
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binary systems, we may constrain the Xc{^'^C) and ^^C(a, 7)^^0 rate more precisely, 
though it requires many calculations for various possible binary pairs. 

One of our purposes for developing the efficient code is to tackle on the evolutions 
just below and above the Fe-core forming SNe. As mentioned in §3.1, such stars 
experience violent shell flashes near the end of their evolution. It is not clear yet 
if such events cause mass loss and shock interactions to be observable. In order to 
follow these stages precisely, it is important for a code to include the acceleration 
term. As mentioned in §2, we will include the term into the YU code in near future. 

The stars just below the critical mass for the Fe-core formation is very interesting 
for another reason. Such a star may become an electron capture supernova (ECSN). 
We would like to construct another progenitor models of ECSNe than the one by 
Ref. [72j) to see if the successful explosion by Ref. 1^1) is model independent or not. 
Since these calculations require high numerical resolution in both time and space,'^ 
we need an efficient code to perform computation. 

We have been working also to develop a new code including the rotation effects 
as in Refs. [T05]) .[106 ]) . [TOTll . [rOS]) .[109 ]) . fTTOl) . This code is based on the YU code and 
aiming for efficient computation as well. Although rotating stellar models in the 
ID formalism have been already calculated by these authors, angular momentum 
transfer is still quite uncertain. Therefore, the evolution of a rotating star is still far 
from complete understanding. Since rotation is inevitable for constructing realistic 
progenitor models for hypernovae and GRBs, and possibly even for normal CCSNe, 
we plan to calculate rotating progenitor models as well. 

The uncertainties in mass-loss rate is the another reason why one needs to com- 
pute several cases to find a better set. Since the knowledge for the rate is still very 
limited, one needs to constrain the rate by using various information including the 
properties of massive stars, supernovae, compact remnants, ISM abundances and 
so on. Unfortunately this may not be easy because other uncertain factors may 
be involved such as rotation, binarity and magnetic field effects. Supernovae may 
provide a key to understand the rate. For example, as discussed in Ref. 20) if one 
can identify SN 2007bi-like event to a CCSN or PISN, we can severely constrain the 
mass loss rates. 

7.2. Nucleosynthesis 

As for our nucleosynthetic works, in §4 we described the success and limitations 
of the simple 'instant energy injection' model. This model reasonably reproduce 
supernova yields up to Zn. Although spherically symmetric hypernova models over- 
produce Fe to lighter elements ratio, such as Fe/0, this problem may be avoided if 
one consider jet-like explosions in 2D for hypernovae even under the instant energy 
injection model. This suggests that the nucleosynthesis up to Zn does not much 
depend on the detail of the central engine. 

On the other hand, nucleosynthesis of r- and weak r-process elements depends 
on the detail of the explosion model. This in turn implies that successful supernova 
and/or hypernova and/or GRB models must produce these elements. 

As described in Ref. HOj) there are observational indications that normal CCSNe 
produce lighter weak r-process elements, Sr, Y and Zr. These elements can be 
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produced and ejected from the hot bubble of normal SNe. However, to produce 
heavier weak r-process elements, Mo and Ru, much larger entropy and neutron rich 
environment is required .'^Sl' It is currently not certain if such an environment can be 
realized in a CCSN. This is because long term simulations of a CCSN showed that 
high-entropy matter ejected from a CCSN is always almost neutral or proton-rich.ED 
As mentioned in §4, it is not clear if all Sr-Zr rich EMP stars are also Mo-Ru rich. 
If this is the case, a CCSN model has to produce Mo-Ru somehow. One possibility 
is the z^p-process fRefs. [TTT]) . [TT2]) . [TT3]) . [TTil) ) not considered in RefJ56l). With this 
process, such weak r-process elements may be synthesized in proton-rich matter 
if entropy is sufficiently high. Thus it is critically important to observationally 
clarify if normal SNe produce Mo-Ru or not. We should note that the results in 
Ref. iGlj) are not obtained by the first principle calculations, but the explosion is 
driven by artificially enhancing the effective neutrino luminosity. The main reason 
that the matter becomes proton-rich is because of the interaction between matter 
and neutrino. If the explosion is driven by the assist of something else, such as the 
rotation energy, neutron rich matter may be ejected. 

The arguments above also applies to the r-process elements. Currently there is 
no clear evidence in the EMP stars that r-process elements are produced in normal 
CCSNe and hypernovae. For example, Zn-rich EMP stars, that we consider them 
polluted by a hypernova, show no enhancement of r-process elements. From the 
results in Ref. [61]) . it seems hopeless to produce r-process elements in normal CCSNe, 
though unknown massive stars should have polluted metal poor r-rich stars, such as 
CS22892-052. Some authors consider NS-NS or NS-blackhole merger systems for 
r-process sites (e.g., Refs . [TT5]) . [TTGI) . ITT7)) . [TTg]) . \TT9\\ . [T20)l ) though it is not clear yet 
if such sites are sufficient to explain the whole r-process abundance in the universe. 
Therefore it is still interesting to explore various SN explosion models to examine if 
they can produce r-process elements or not. 
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